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^ [Abstract 

The incorporation of particle inertia into the usual mean field theory for particle aggregation and fragmentation in 
fluid flows is still an unsolved problem. We therefore suggest an alternative approach that is based on the dynamics of 



^ individual inertial particles and apply this to study steady state particle size distributions in a 3-d synthetic turbulent flow. 
^ We show how a fractal-like structure, typical of aggregates in natural systems, can be incorporated in an approximate 
•way into the aggregation and fragmentation model by introducing effective densities and radii. Based on this model 
■we investigate numerically the impact of three different modes of fragmentation: large-scale splitting, where fragments 
have similar sizes, erosion, where one of the fragments is much smaller than the other and uniform fragmentation, where 
— I'all sizes of fragments occur with the same probability. We flnd that the steady state particle size distribution depends 
^ strongly on the mode of fragmentation. The resulting size distribution for large-scale fragmentation is exponential. As 
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some aggregate distributions found in published measurements share this latter characteristic, this may indicate that 
large-scale fragmentation is the primary mode of fragmentation in these cases. 

'Keywords: inertial particles, aggregation, fragmentation, fractal aggregates 
PACS: 05.45.-a, 47.52. -hj, 47.53.+n 



1. Introduction 

In recent years there has been a great effort to in- 
vestigate the advection of inertial particles in fluid flows 
0, S S S 0] • Understanding the behavior of inertial 
particles like aggregates, dust or bubbles moving in in- 
tompressible flows plays an important role in such diverse 



fields as cloud physics 
sediment dynamics 



X 



engineering [8] , marine snow and 

^ as well as wastewater treatment 

The dynamics of these particles is dissipative, which 
leads on the one hand in most flows to a preferential ac- 
cumulation in certain regions in space, i.e. on attractors, 
and on the other hand to a behavior that is very different 
from passive tracers 12 , 13|, [l3| . Previous studies concen- 
'trated mainly on non-interacting particles, in spite of the 
fact that accumulation leads unavoidable to mutual inter- 
actions of different kinds. 

It is well known that as a result of collisions between par- 
ticles, aggregates can be formed that consist of a large 
number of primary particles. In many areas of science 
the formation of such aggregates and their break-up due 
to shear forces in the flow plays a very important role, 
e.g. in sedimentation of particles in oceans and lakes [l5l |. 
chemical engineering systems such as solid-liquid separa- 
tion 0, 17 [, aggregation of marine aggregates ^ISi] and 



flocculation of cells [19| . 
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Most approaches of aggregation and fragmentation mod- 
els are based on the pioneering work of Smoluchowski [20j 
and use usually a mean field approach with kinetic rate 
equations to model these processes (see e.g. Jackson Q). 
However, for particles with inertia a field theory for the 
particle velocity has not yet been formulated. The exis- 
tence of caustics, meaning that the dynamics of inertial 
particles would lead at some points to a multi- valued par- 
ticle velocity field [21, 22|, has prevented such an approach 
so far. While attempts have been made to incorporate par- 
ticle inertia in approximate ways in a mean field approach 
23 . 2^ . no completely satisfying solution has been found 



yet. 

Here we therefore choose a different, individual particle 
based approach, where the dynamics of finite size parti- 
cles are taken directly into account. The approach has 
recently been introduced in Ref. [l^, and discussed in 
more detail with respect to different fiows in [2^. In the 
present study we adopt this approach to study the long- 
term behavior of particle size distributions that develop 
from a balance between aggregation and fragmentation. 
In particular, we examine the influence of fragmentation 
and aggregate structure on these size distributions. 
In most previous works the particles were considered to 
be spheres with a specific density. In many realistic cases, 
for example for marine aggregates, this is a crude approx- 
imation. The complex structure of particles can have a 
great influence on particle dynamics as well as aggregation 
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and fragmentation processes. Both the actual motion of 
aggregates and the probabilities for aggregation and frag- 
mentation are influenced by the structure of the particles. 
In the context of a mean field approach, a complex par- 
ticle structure has been incorporated in the past in terms 
of a density modification for the particles, e.g. by Kra- 
nenburg 21] or Maggi et al. [2^. However, so far there 
are very few attempts to treat this problem for inertial 
particles in a fiow. Wilkinson et al. [29'| used a model for 
fractal particles in an aggregation model for dust particles 
during planet formation. Our present work expands the 
consideration of spherical particles to model more realistic 
aggregates. We focus specifically on the problem of aggre- 
gation and fragmentation of marine aggregates, who can 
often be described as having a fractal-like structure 27 1. 
By this we mean that on average there exists a power- 
law relationship between the characteristic length and the 
mass of such aggregates. The exponent of the power-law 
is called the fractal dimension. Such a characterization in 
terms of a fractal dimension leads to a modification of the 
radii and effective densities of the aggregates compared to 
a solid sphere of the same mass. Nevertheless, we still 
treat them as effectively spherical for the particle motion, 
allowing us to apply the Maxey-Riley equations of motion 
(30I with modified parameters. 

In this work we choose a parameterization of our model 
for the case of a suspension of marine aggregates in the 
ocean. In this way we can study our modeling approach 
for a specific case, but we emphasize that our model is a 
general one that can in principle be used for a wide range of 
applications where aggregation and fragmentation of solid 
particles plays a role. 

Since the fractal dimension of marine aggregates can vary 
greatly in natural systems [28|, we examine its effect on 
the steady state particle size distributions in our model. 
We find that while the shape of the size distributions does 
not depend strongly on the fractal dimension, the average 
particle size and relaxation time towards the steady state 
depend strongly on this parameter. 

Even though to a certain extent methods from dynamical 
systems theory can usefully be applied, we mention that 
the entire problem is much more complex than that of any 
usual dynamical system. While particles of a single size 
move on specific attractors, aggregation and fragmentation 
lead to repeated transitions from one attractor to another 
one, depending on the aggregate size. The skeleton of the 
new dynamics is therefore a superposition of the different 
attractors, with transient motion in between. The struc- 
ture of the individual attractors and their superposition in 
turn influence the aggregation probabilities due to differ- 
ent local concentrations of particles. Fragmentation is also 
affected by the particle dynamics, because shear forces can 
be locally different in the flow. Therefore, break-up may 
depend on whether an attractor for a certain particle size 
lies in a region with high shear or not. 
We show that the combination of aggregation and frag- 
mentation of fractal-like aggregates, superimposed on in- 



ertial advection dynamics, leads to a convergence to a 
steady state in the particle size distribution. This steady 
state is unique for a given set of parameters. Mainly, we 
compare three different types of splitting, uniform frag- 
mentation, erosion and large-scale fragmentation. These 
splitting modes differ in the size of the fragments that 
are created during break-up. While erosion leads to one 
large and one relatively small fragment, large-scale frag- 
mentation leads to two fragments of similar size. We find 
that the transient dynamics as well as the size distribu- 
tion in the steady state depend strongly on the splitting 
mode. The steady state size distribution found for large- 
scale fragmentation conforms best to observation reported 
in the literature for the break-up of marine aggregates in 
tidal areas [sij, indicating that this may be the primary 
mode of fragmentation in these cases. 
Section [2] describes the complete model for advection, ag- 
gregation and fragmentation that is used in this work. The 
equations of motion for heavy spherical particles (Stokes 
equation) are used, but with modified parameters to take a 
fractal-like structure into account. Rules governing the ag- 
gregation and fragmentation are introduced. Finally, the 
model is applied to a simple 3-d synthetic turbulent fiow 
field. 

Section [3] then presents a complete analysis of the in- 
fluence of all major system parameters on the resulting 
steady state size distributions, the average aggregate size 
in steady state and the relaxation time towards the steady 
state. Namely, these parameters are aggregate strength, 
fractal dimension of the aggregates and total particle mass 
in the system. 

Section 3] contains a discussion of the limitations of the 
model and the conclusions. 



2. Advection, Aggregation and Fragmentation 
Model 

In this section we will present the modeling approach 
used in this study, that describes the motion, aggrega- 
tion and fragmentation of flnite-size particles. Firstly, the 
equations of motion used for the advection of spherical 
particles heavier than the surrounding fluid are presented. 
Secondly, a model to account for the fractal-like structure 
of real aggregates is described. Thirdly, a full model to 
include aggregation and fragmentation in this context is 
introduced. Finally, a simple 3-d synthetic turbulent flow 
fleld is chosen, that will be used to study the aggregation 
and fragmentation model in detail. 

2.1. Equations of Motion for Spherical Particles 

For simplicity, we consider all primary (smallest, un- 
breakable) particles to be spherical and denser than the 
ambient fluid. We emphasize that the equations of motion 
presented here for spherical particles will in the following 
also be used to describe the motion of ag greg ates which 
usually can not be assumed to be spherical [21]. However, 
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to account for properties related to the fractal-like struc- 
ture of aggregates some modifications to the equations of 
motion (in form of modified parameters) will be introduced 
in the next section. 

Finite-size particles usually do not follow exactly the mo- 
tion of the surrounding fluid, instead inertia effects lead 
to deviations of the particle motion from that of the fluid. 
For small Reynolds numbers the equations of motion for 
spherical particles of finite size are the Maxey-Riley equa- 
tions [s^. While inertia effects can be fairly small for 
the primary particles in the case of marine aggregates (see 
section [2761) . the influence of particle inertia increases with 
aggregate size and can become quite important for larger 
marine aggregates. 

In reality every particle produces perturbations in the flow 
that decay inversely proportional to the distance from the 
particle [321 • In this work we keep the particle concentra- 
tion n low enough to be in a diluted regime. For particles 
of radius r and density pp moving coherently within the 
dissipative scale Ip of a flow the feedback from the parti- 
cle motion on the flow can be neglected if nrt^ <^ 1 [33l |. 
Particle-particle interactions mediated by flow perturba- 
tions are neglected, see the discussion in section ITT] 
Assuming that the difference between the particle velocity 
V(i) and the fluid velocity u — u(X.{t),t) at the position 
X(t) ^ {Xi{t),X2{t),X3{t)) of the particle is sufficiently 
small, the drag force is proportional to this difference. This 
is called Stokes drag. With these restrictions the dimen- 
sionless form of the governing equation for the path X(t) 
of such a particle under the influence of drag and gravity 
can then be approximated from the Maxey-Riley equations 
as: 

V=^(u-V-Wn), (1) 

where n is the unit vector pointing upwards in the vertical 
direction (which is the X2-axis in this study). 
The particle Stokes number St, i.e. the ratio between par- 
ticle response time and flow time scale is defined as 

St = {pp2r')/{9fiTf) (2) 

and the dimensionless settling velocity in a medium at rest 
W is defined as 

W^2r^pp^PF)Tfg/{9t,lf) . (3) 

Here, pp and p are the fluids' density and dynamic viscos- 
ity and I f and Tf are characteristic length and time scales 
of the flow. 

2.2. Fractal-like aggregates 

When looking at real aggregates they are typically not 
solid spherical particles but instead can have a complex 
structure as they consist of a number of primary parti- 
cles. In this model the primary particles are assumed to 
be solid spherical particles, following the equations of mo- 
tion as described in the last subsection. All aggregates are 



assumed to consist of an integer number of such spheri- 
cal primary particles. The description of the motion of 
an aggregate with a complex structure is still an unsolved 
problem. Therefore we only consider the influence of the 
structure of the aggregates on their size and effective den- 
sity. Here, we assume that aggregates have a fractal-like 
structure, meaning that there exists a power-law relation- 
ship between the characteristic length and the mass of such 
aggregates. The structure of the aggregates can then be 
characterized by a fractal dimension c?/ < 3 Their 
size can still be defined approximately by a radius, that 
can be considered as the characteristic length scale of the 
aggregate. This radius of an aggregate that consists of 
a primary particles and has a given fractal dimension d/ is 
derived in the following. We emphasize that the number 
of primary particles a in an aggregate is here also used 
as an index to describe a quantity, e.g. the radius or the 
volume, that corresponds to an aggregate consisting of a 
primary particles. 

The solid volume, i.e. the volume of an aggregate that is 
filled with particulate matter follows from the definition of 
the fractal dimension df (see e.g. ^) as 

solid — C(lfrc/ , (4) 

where c^f is a proportionality constant that can depend 
on df. As mass conservation must be fulfilled we get 

Va.soiid = aVi = a^TTrl , (5) 

where ri and Vi are the radius and volume of a primary 
particle, respectively. The proportionality constant Cdf 
can be derived from Eqs. (|H) and (O by setting a — 1 

M 

Cdf = gTrr^ . (6) 

In combination with Eqs. (|1]) and ([S]) this leads to 

r„ = a^^^'fa, (7) 

for the radius of an aggregate. It is evident that due to 
the fractal-like structure the radius is greater than for 
a completely solid particle of the same mass. 
A part of a fractal-like aggregate, i.e. of the total volume 
encased by an aggregate 

Va, total = ^Trr^ , (8) 

is not filled with matter but with the surrounding fluid. 
The aggregate density therefore decreases with an increas- 
ing number a of primary particles in the aggregate. From 
mass conservation it follows 

Pa V'a, total — Pi V'a, solid + PF(Va, total " solid) (9) 

Solving Eq. (jO]) for pa and substituting Vq, total and T4,_soiid 
leads to 

Pa = PF + (pi - PF)a'"'/''^ . (10) 
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Going back to the equations of motion we now see that 
as far as the particle dynamics are concerned, a first ap- 
proximation for the fractal-hke aggregates is to treat them 
as spheres with an increased radius compared to sohd 
spheres of the same mass, but a reduced density pa, be- 
cause of the fluid encased in their spherical volume. 
This leads to a modification of the particle Stokes number 
St and dimensionless settling velocity W for fractal- like 
aggregates in the equations of motion ([T]) , when compared 
to a solid sphere: 



W„ = 2Tl{p^-pF)Tfg/{%tilf) 



(11) 
(12) 



For fractal- like aggregates these parameters replace St and 
W in Eq. ([T]), leading to a motion with different param- 
eters for aggregates with different numbers a of primary 
particles. 

2.3. Aggregation model 

The physical, chemical or biological process that leads 
to aggregation of particles can vary from case to case and 
is not examined in detail here, as this is beyond the scope 
of this study. Instead, a general model is used. The only 
assumption is that during a collision particles can some- 
how stick together and form an aggregate. No detailed 
mechanism of sticking is considered. 

Whenever the distance between two particle centers be- 
comes smaller than the sum of their radii, these two par- 
ticles aggregate immediately, creating a new particle that 
replaces the two old particles. For two particles of radius 
Oq. and Octj the new number of primary particles after ag- 
gregation is obviously Ofncw = Q^i + ctj, leading to a new 
radius aa„„„, which can be derived from Eq. ([7]). The 
position of the new particle is the center of gravity of the 
old particles. The velocity of the new particle follows from 
momentum conservation. 

2.4^. Fragmentation model 

The physical process leading to fragmentation can 
vary as well. While there are detailed models for the 
fragmentation of e.g. water droplets .37], the mecha- 
nism of fragmentation of marine aggregates is not well 
understood. Even experiments give no clear indication 
how fragmentation of such aggregates occurs in detail. 
Therefore only a few theoretical approaches exist for this 
process 38, 39]. In this work a model for fragmentation 
is used that is only based on some very general properties 
of the aggregates involved. 

In the following the fragmentation is described in two 
parts, that need to be clearly distinguished. Firstly, a 
splitting condition^ that determines if a fragmentation 
event takes place and secondly, a splitting rule, that 
determines how fragmentation takes place, are defined. 



2.4.1. Splitting condition 

The splitting condition describes if fragmentation of an 
aggregate takes place. Generally, only particles that are 
composed of more than one primary particle can fragment. 
The break-up of an aggregate occurs when the hydrody- 
namical forces i^hyd acting on the aggregate exceed the 
forces Fagg holding the particles in the aggregate together. 
The criterion for breakup can therefore be expressed as 



^hyd/^kgg > const. 



(13) 



For aggregates with a fractal-like structure, consisting of a 
number of solid spheres, the hydrodynamical forces in the 
dissipative range where viscous forces dominate is propor- 
tional to the shear force integrated over the surface of the 
aggregate ^0]. For a fracture at a distance C-r^, C S [0, 1[ 
from the equator of the aggregate this results in 



(14) 



where S is the shear rate in the flow. 
For a porous aggregate Ruiz et al. [4lj give the force 
Fagg holding an aggregate together as proportional to the 
area of constituent matter in a cross-section of the ag- 
gregate. Ruiz et al. then related Ki, solid to the porosity 
of an aggregate. Here, we instead rewrite this relation- 
ship in terms of the fractal dimension df. For fractures 
across the equator of an aggregate the area is proportional 



to V 



2/3 

, solid 



oc r 



2dj/3 



For a fracture at a distance 



C-ra, C G [O; 1] the area of constituent matter in the cross- 
section is reduced. It equals the area of a cross-section 
across the equator of an aggregate with decreased radius 
(1 — C^y^^ra. In general we get 



Fagg a (1 - Cr^^'rl {rjr.r^'^' 
The splitting condition ([T^ then becomes 



S- 



(1 - er/"" fr_^ 

(1-C) A^i 



> 1, 



(15) 



(16) 



where the proportionality constant 7 is determined by the 
force holding the primary particles in an aggregate to- 
gether. It is therefore a measure of the aggregate strength. 
Solving for the shear rate S and using Eq. ([7|) leads to 



= (i^c^y^ 2/3-2/d, 
^ (1-0 



(17) 



It can be seen that for a fractal dimension d/ < 3 the 
critical shear force required to break up an aggregate de- 
creases with the aggregate size, i.e. larger aggregates are 
less stable than smaller ones. Additionally, the critical 
shear required for fragmentation is smallest for fractures 
across the equator of an aggregate [C, = 0) and increases 
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with increasing distance from the equator. 
The shear force S is given by 



(1/2) 



— I 2 ^ Sij Si- 



(18) 



where Sij = \ + ^y") is the rate-of-strain tensor in 

the flow. 



a) Large-scale fragmentation 



2.4.2. Splitting rules 

The sphtting rules describe how an aggregate wih spht, 
when the sphtting condition is met. During fragmenta- 
tion only particles whose mass is an integer multiple of the 
mass of a primary particle are created. This means that, 
even though they have become part of some larger aggre- 
gates, primary particles can never be broken up. Only the 
connection among each other can break. Different size dis- 
tributions of the fragments are possible. 
When a splitting condition is met, an aggregate consisting 
of aoid primary particles is split into 2 fragments with the 
number au of primary particles of each fragment being a 
random fraction of the original number aoid- 
Typically, one distinguishes between two different mecha- 
nisms of fragmentation 42| . For each mechanism fragmen- 
tation occurs on average at a different distance from the 
equator of the aggregate, leading to different distributions 
of the fragments a^. This can be expressed in different 
values for ^, the fraction of the distance from the equator 
of the aggregate where fragmentation is assumed to take 
place (see Sec. I2.4.ip . Large-scale fragmentation happens 
when an aggregate is 'pulled apart' somewhere close to the 
equator, leading to fragments of similar size. This is char- 
acterized by = 0. Erosion happens when shear forces 
act closer to the edge of an aggregate [43[, implying that 
< C < 1 . In this case only few primary particles are split 
off from the aggregate (see Fig. [T]). 

Here we will compare three different fragmentation 
modes, with different distributions for the number of pri- 
mary particles in the fragments. First a large-scale split- 
ting rule, second an erosion splitting rule, and third a uni- 
form splitting rule. In most realistic cases one expects 
different fragmentation modes to appear together, even 
though with slightly different probabilities. However, here 
we apply these fragmentation modes separately to deter- 
mine their individual influence on the aggregate size dis- 
tribution. 

1. For both large-scale fragmentation and erosion we as- 
sume that there is a preferred distance Qra, C ^ 1[ 
from the center of the original aggregate where frag- 
mentation occurs. The two mechanisms are then 
characterized by different values of Q. A fracture 
at a distance C,ra leads to a fragment containing 
a fraction Vpr/Va^totai of the original volume, where 
Vfj. — 7rrj^(| — C + Since C, is assumed to be only 
the average distance of a fracture from the equator of 



b) Erosion 



Figure 1: Sketch of different splitting meclianisms: a) large-scale 
fragmentation where an aggregate is pulled apart along the equator, 
b) erosion where small parts are split off from the aggregate surface. 



the aggregate, for each fragmentation event we choose 
the number of primary particles in the fragment ran- 
domly from a Gaussian distribution, centered around 
Q!oid • Vh /Va^totai- This allows for a certain variation 
of the fragment size, meaning that 



ai 



3 I3 ^ 3 



C "old 



(19) 



where ^ is a random number from a normal distri- 
bution with mean and standard deviation ap and 
the right-hand side of equation (TT^ is rounded to- 
wards the nearest integer. We note that our results 
do not depend strongly on the choice oi ap, here we 
choose Up = 0.2 which results in a typical variation of 
one primary particle between fragments coming from 
identical aggregates. As an additional restriction it is 
required that 1 < ai < aoid, otherwise a new random 
number ^ is chosen. 

For large-scale fragmentation, aggregates are assumed 
to break along the equator into two fragments of simi- 
lar size, which corresponds to 0. For erosion, frag- 
mentation is assumed to occur at a distance from the 
center, leading to one aggregate being much smaller 
than the other. This corresponds to < C < 1. 
Here, we choose C = 0.6 for erosion, which leads to 
smaller fragments containing on average 10 percent 
of the mass of the original aggregate. Similar results 
for erosion were found for other values of C- However, 
if ^ becomes too large particles will no longer frag- 
ment because the critical shear required to break off 
a fragment increases greatly as C — > 1. 
2. In the uniform splitting rule the number of primary 
particles for the first fragment is chosen from a uni- 
form distribution in the interval / — [l,aoid[- The 
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uniform splitting rule is a simplified model for the full 
case where both large-scale fragmentation and ero- 
sion of an aggregate can happen. However, all sizes 
of fragments occur here with the same probability. 

For all three cases the remaining aggregate consists of 
a2 = ctoid — cti primary particles. 

Whenever a particle is split according to one of the split- 
ting rules, all parts keep the velocity of the original parti- 
cle. That way momentum is conserved. The first fragment 
remains at the position X = {Xi, X2, X-s) of the original 
particle. The center of the other fragment is placed along 
a line segment in a random direction, so that for the two 
fragments the distance equals the sum of the radii. 
For each fragment, the splitting condition is checked again 
and if it is met, the whole process is repeated until no 
fragment fulfills the fragmentation condition. This leads 
to a splitting cascade and aggregates can break up into 
more than two fragments, if the aggregate is large enough 
or shear forces are strong enough. This is consistent with 
experimental observations of marine aggregates that larger 
particles tend to break into more fragments 44]. In that 
way ternary, quarternary and other splitting types besides 
binary splitting naturally appear in this model. 
Here, large-scale fragmentation and erosion are treated as 
two separate processes to study the influence of fragmen- 
tation at certain distances on the aggregate size distribu- 
tions in the steady state. In reality, aggregates will break 
with certain probabilities at certain distances C from the 
center but there will be no two separate processes. There- 
fore, depending on the probability distribution for frag- 
mentation at a certain distance one can expect different 
combinations of the steady state size distributions found 
in this work. The uniform fragmentation rule represents 
one such possible combination, where fragmentation at all 
distances C appears with the same probability. 

2.5. Fluid Flow 

As a fluid velocity field we consider synthetic turbulence 
in the form of a space-periodic, isotropic and homogeneous 
Gaussian random flow [i^, since it allows us to perform 
long-term simulations at reasonable computational costs. 
We use a smooth, incompressible flow since we focus on 
effects typically taking place on scales smaller than the 
Kolmogorov scale of a turbulent flow. 

The flow is written as a Fourier series 

u{X,t)= Yj Hk,t)e'^''-^ , (20) 

fcGZ<i\{0} 

where u{k,t) G are the Fourier components, with the 
property u{—k, t) — u*{k, t) because u{X, t) is real- valued. 
The star denotes complex conjugation. By taking for 
u{k,t) the projection of a different vector v{k,t) G 
onto the plane perpendicular to the wave vector k, incom- 
pressibility is ensured. The vector v{k,t) is assumed to 



be an Ornstein-Uhlenbeck process. It is a solution of the 
complex-valued stochastic differential equation 

S = -^^{k)vdt (j{k)dW , (21) 

with £^{k),<T{k) gM, where dW is a d dimensional complex 
Wiener increment. The parameters ^(k),a{k) need to be 
chosen in such a way that the flow u(x, t) reproduces some 
features of a real turbulent flow, in this case we choose the 
energy spectrum in the dissipative range of a turbulent 
flow. Here we use the exponential spectrum suggested by 
Kraichnan 

E{k) = C ■ {2Trklf/Lf exp(-/3[27rfc///L]) , (22) 

with /3 = 5.2 [4^ and a suitably chosen normalization 
constant C. The constant / / is the length scale of coherent 
structures in the flow and L is the spatial period of the 
flow. We choose ^(fc) = 1/t/ and cr(fc) = ^E{k)/Tf. The 
constant t/ is then the correlation time of the flow. The 
normalization constant is chosen so that a desired value 
of the mean shear rate {S) is obtained. The flow is then 
characterized by the correlation time Tf, the correlation 
length I f and the mean shear rate {S) . 

If a fluid velocity field with few Fourier modes is chosen, 
no interpolation of the velocity at particle position is re- 
quired, since it can be calculated from direct summation 
of the Fourier series. This allows for a resolution of the 
fine structures of the particle distribution in space. 

2. 6. Implementation 

Next, we will describe some specifics of the numerical 
implementation and the system parameters used in this 
work. 

For the simulations in this work we choose particle proper- 
ties similar to those of marine aggregates in coastal waters. 
The primary particles considered in this model have a ra- 
dius ri — 4/im, density pi — 2.5 x lO^kg/m^ and mass 
= pi|7rrf. The relevant characteristic length scales 
for marine aggregates in coastal areas of the ocean are 
typically the Kolmogorov scales. Shear rates in coastal 
areas can be of order ls~^, leading to Kolmogorov length 
and time scales oi If = 1mm and Tf = Is, respectively 
[27I I. Using these scales to make the equations of motion 
of the particles dimensionless leads to a Stokes parame- 
ter of Sti = 10^^ and a dimensionless settling velocity of 
Wi — 0.1 for the primary particles. 

The aggregate strength parameter 7 is fixed at 7 = 8, un- 
less otherwise mentioned. 

The number of aggregates N{t) changes over time due to 
aggregation and fragmentation leading to a distribution of 
aggregates of different radii in the flow. However, the to- 
tal mass M — "Y^^J^i aivni remains constant during one 
simulation. As initial condition we take 10^ primary par- 
ticles and no larger particles. Furthermore particles are 
uniformly distributed over one periodic cell of size of 
the conflguration space, with velocities matching that of 



6 



the fluid. This choice fixes the total mass of the system to 
be M = lO^mi. For the flow we choose a periodic cube 
with L = Alf, so that we obtain a volume fraction of about 
0.4 X 10"^. 

The fractal dimension df oi marine aggregates varies be- 
tween approximately 1.5 for very open, fragile aggregates 
like marine snow in the open ocean and approximately 2.5 
for stronger, compact aggregates. The average is typically 
around 1.9 — 2.0 (see e.g. [la]), therefore in the following 
we choose df = 2.0 unless otherwise mentioned. 
As a flrst approximation the aggregation and fragmenta- 
tion processes are assumed to have no effect other than 
to change the size of the particles and the effective den- 
sity, and therefore do not directly influence the motion of 
the particles. Hence all three aspects, motion, aggrega- 
tion and fragmentation that define the whole system can 
be modeled separately. Aggregation is checked constantly 
during the integration, whereas fragmentation is applied 
after every time step of the system. 

1. All particles move in the flow for one time step dt 
using the equations of motion described in Sec. 12.11 
We emphasize at this point that each aggregate size 
is characterized by different values of St a and Wa , so 
that the motion of aggregates of different size is gov- 
erned by the same equations but with different pa- 
rameters. 

The length of the time step dt needs to be chosen 
small enough so that the simulation result becomes 
independent of this values, here we found dt — 0.01s 
to be sufficiently small. 

Because of the spatial periodicity of the flow, all par- 
ticle dynamics will be folded back onto one cell in 
the flow, using periodic boundary conditions. Usu- 
ally particles do not stay in one cell of the flow, i.e. 
they are not suspended in the flow. Instead parti- 
cles will generally fall downwards through the flow, 
if they are heavier than the fluid [47]. This means 
that folding the dynamics of the particles back onto a 
single cell is only a convenient way to visualize an in- 
finitely extended system and does not completely mir- 
ror what one would see in a comparable experiment. 
However, if particles are initially distributed homoge- 
neously over the whole configuration space, the total 
particle mass in each periodic cell remains the same 
over time. Therefore even if aggregation and frag- 
mentation are included, it is sufficient to restrict our 
studies to one unit cell with periodic boundaries. 

2. Particles aggregate upon collision, i.e. if their dis- 
tance becomes equal to the sum of their radii. To 
ensure that no collisions are missed, we use an effi- 
cient event-driven algorithm for particle laden fiows 
(cf. Sigurgeirsson et al. for details). Computa- 
tionally, the aggregation process is the most costly 
component of the simulation. In particular, the naive 
approach to check which particles are colliding in- 
volves looping over all pairs of particles and therefore 



scales as O(iV^), where N is the number of particles. 
Therefore, here a linked-list algorithm [4^, sometimes 
also called link-cell algorithm is used to compute the 
distance between particles. The configuration space 
is divided into grid cells of size e, where each grid cell 
stores information on which particles it contains. The 
looping over particle pairs to calculate their distance 
is then done only over particles in a given grid cell and 
the neighboring cells. If the grid cell size e is small 
enough (but larger than the largest appearing parti- 
cle size) the link-cell algorithm scales as 0{N) and is 
thus much faster than the naive approach. 
3. After each time step dt particles can fragment if the 
shear at their position exceeds a critical value. If that 
is the case, new fragments are created according to 
the rules described in section [2.4.21 

We note that at first glance it looks like aggregation 
and fragmentation are treated very different, in particular 
aggregation seems to be independent of the aggregate 
strength 7 in this model. However, this is not the case. 
Initially all particles that come into contact aggregate, i.e. 
here the aggregation probability upon collision is equal to 
one. But when looking at aggregation and fragmentation 
together over one time step dt it is in fact smaller than 
one because some aggregates that just formed during this 
time step will break up again. These are the aggregates 
where the aggregate strength 7 is not strong enough to 
hold the aggregate together. This means after one time 
step only some of the particles that came into contact will 
actually have aggregated and this number will depend 
on the aggregate strength 7. This means that both 
aggregation and fragmentation probabilities depend on 
the same aggregate property, which one would expect in 
reality. 



3. Simulation results 

In the following section we will present simulation 
results using the model described above, to determine the 
infiuence of the different splitting rules on the resulting 
particle size distributions. As the parameters used in 
the model system can vary greatly in natural systems, 
we examine the sensitivity of the system with regard to 
the following parameters: aggregate strength 7, fractal 
dimension df and total number of primary particle N. 

3.1. Measured quantities 

From previous works it is known that the balance of 
aggregation and fragmentation typically leads to a steady 
state of the particle size distribution [H3|. This follows 
from the fact that normally aggregation dominates for 
small sizes, whereas fragmentation is the dominant pro- 
cess for large sizes. In addition to studying this size dis- 
tribution of the particles in the steady state, we introduce 
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different measures to characterize first the approach to the 
steady state and then the steady state itself. 
To follow the convergence of the system towards a steady 
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Figure 2: Average number of primary particles as a function of 
time for aggregate strength 7 = 9 and aggregate fractal dimension 
df = 2.0 for (a) large-scale fragmentation (b) erosion and (c) uniform 
fragmentation. 

state, we use two different quantities. The first quan- 
tity that we measure during the simulations is the average 
number of primary particles per aggregate, defined as 



(23) 



In the context of our model 2 {a{t)) corresponds to the 
'mean equivalent circular diameter' that is often used as a 
measure in experiments with marine aggregates |3lj. We 
will use this quantity as a first estimate whether the par- 
ticle size distribution has converged to a steady state and 
to follow the evolution of the particle size distribution to- 
wards the steady state. 

A second quantity of the aggregation and fragmentation 
process that may be experimentally measured is the time 
it takes to reach the steady state. Especially in technical 
applications this can be an important quantity, where pro- 
cesses need to be timed appropriately to allow for a smooth 
work flow. Here we introduce a measure for this relaxation 
time in our model and show how different system param- 
eters influence this time to reach the steady state. 
We define the relaxation time Tco as 



dt 



(24) 



{a{t)) is a running (time)-average of the average number 
of primary particles in an aggregate. It is used to remove 



fluctuations due to the periodic changes in the flow. This 
definition of the relaxation time is analogous to the defini- 
tion of the correlation time for stochastic processes as the 
integral over the autocorrelation function. For example, 
for an exponential relaxation process oc e^*/*" Eq. ([M)) 
leads to the expected result of Too — tR- However, Eq. 
(|24p stays an appropriate measure for more irregular re- 
laxation processes. 

As a simple measure to characterize the steady state of 
the particle size distribution we use the average number of 
primary particles in an aggregate in the steady state that 
is defined as 

aoo = lim {a{t)) . (25) 

3.2. Approach to a steady state 

First, we use the average number of primary particles 
per aggregate to follow the convergence of the systems 
towards a steady state for all three splitting rules (Fig. 
[2]). Our initial condition is always a uniform distribution 
of primary particles. 

Initially, aggregation leads to a fast increase in the average 
number of primary particles per aggregate similar for all 
splitting rules. Then fragmentation sets in and a balance 



c lO" 



SlO" 



lO"- 



a 10"- 



1 1 


large-scale — 
uniform 


-o - - 


■ \ 


erosion — 


-T 








■ \ "''n 






\ ^ 






1 \ 1 







^ -4 

10^ 

23456789 
# of fragments 

Figure 3: The histogram shows the percentage of fragmentation 
events for the number of fragments created in that event, for the 
same simulation run as shown in Fig. [2] Large-scale fragmentation 
leads typically to the smallest number of fragments, while erosion 
typically generates most fragments. 

between aggregation and fragmentation is reached, with 
a different steady state average particle size for the 
different splitting rules. fluctuates over time, due 

to the statistical fluctuations in the flow. Large-scale 
splitting leads to the highest average number of primary 
particles per aggregate, erosion to the lowest and uniform 
fragmentation is in between. This can be intuitively 
understood, since for erosion typically more fragments 
are created than for large-scale fragmentation (see Fig. 
[3]). When a particle gets eroded, one of the fragments is 
usually close to the same size as the original aggregate. 
This leads to a high probability that this fragment will 
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break again and therefore in many cases fragmentation 
will not be binary, but many fragments will be created. 
For large-scale fragmentation, aggregates will typically 
break only once, since both fragments are much smaller 
than the original aggregate. 

In general, it is less likely in the case of large-scale frag- 
mentation that a large number of fragments is created. 
This leads in the mean to a larger average aggregate size 
than for erosion. 
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Figure 4: Influence of aggregate strength 7. (a) Average number of 
primary particles per aggregate in steady state Ocxd and (b) relaxation 
time Too for the approach to the steady state for different values of 
the aggregate strength 7. Error bars are obtained from an ensemble 
of 5 different realizations of the carrier flow. 

The different splitting rules lead to very different 
distributions (cf. Fig. [S|). Large-scale fragmentation 
creates a distribution with a single peak at intermediate 
radii and no particles of the smallest sizes. The right 
hand side of the aggregate size distribution for large-scale 
fragmentation follows approximately an exponential 
decay. The size distribution found for large-scale frag- 
mentation corresponds well to those observed for marine 
aggregates [sH where exponential size distributions have 
also been reported. This may indicate that large-scale 
fragmentation is indeed the primary mode of break-up for 
many marine aggregates, as proposed in some works (see 
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Figure 5: Histogram of the particle size distribution. Number of 
particles versus the radius a for different values of 7 for (a) large- 
scale splitting (b) erosion splitting and (c) uniform splitting. 

By contrast, the size distribution for erosion has two 
different regimes, with a sharp maximum at the smallest 
aggregate size and a slower decaying tail at larger aggre- 
gate sizes. 

Uniform splitting, where both larger and smaller fragments 
are created leads to a plateau in the size distribution at 
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smaller aggregate sizes and an exponential decay towards 
larger aggregate sizes. 

Many of the system parameters that appear in our model 
can vary much in natural systems, in particular the aggre- 
gate strength, the number of primary particles involved 
and the fractal dimension of aggregates. In the following 
we therefore examine the sensitivity of our results to these 
parameters. 

3.3. Influence of aggregate strength 

To determine the influence of the forces holding the 
aggregates together on the resulting steady state size 
distribution, aoo is computed for different values of 
the aggregate strength 7. aoo increases with 7 for all 
fragmentation rules. The increase is fastest for large-scale 
fragmentation and slowest for erosion, (see Fig. HJa)). 
The relaxation time Too as a function of 7 is shown in Fig. 
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Figure 6: Influence of the total number of primary particles A^o- (a) 
Average number of primary particles per aggregate in steady state 
ctaa and (b) relaxation time r^a for the approach to the steady state 
for different values of the totalnumber of primary particles A^o ■ Error 
bars are obtained from an ensemble of 5 different realizations of the 
carrier flow. 

m^b). For all three fragmentation rules, the relaxation 
time is independent of the value of 7. 
The relaxation time is defined relative to the average 
number of primary particles in steady state aoo and 



Figure 7: Histogram of the particle size distribution. Number of 
particles versus the radius a for different values of the total number of 
primary particles A^o for (a) large-scale splitting (b) erosion splitting 
and (c) uniform splitting. 

When looking at the particle size distribution in steady 
state (Fig. [SJ, the difference between the three fragmenta- 
tion rules is again clearly visible. However, this difference 
does not seem to depend on the value of the aggregate 
strength 7, as the distributions for each fragmentation rule 



10 



remain qualitatively the same for different 7, but getting 
wider with increasing aggregate strength. 

3.4- Influence of the volume fraction 

To determine the influence of the volume fraction, i.e. 
the total number of primary particles iVo in the system, 
on the resulting size distribution, aoc is computed for 
different values of Nq. aoo increases with increasing 
number of primary particles, due to the increased number 
of collisions. Again, the average number of primary 
particles per aggregate in the steady state aoo is largest 
for large-scale fragmentation and smallest for erosion. 
The increase in aoo with A^o is almost linear for all three 
fragmentation rules. 
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Figure 8: Influence of the fractal dimension df of the aggregates, (a) 
Average number of primary particles per aggregate in steady state 
OcxD, the inset shows the relative average size of aggregates in the 

steady state Too I fx = ctxf^^ of the aggregates as a function oi df 
and (b) relaxation time Too for the approach to the steady state for 
different values of the fractal dimension df o{ the aggregates. Error 
bars are obtained from an ensemble of 5 different realizations of the 
carrier flow. 

The relaxation time (see Fig. Elb)) decreases for 
increasing M for all three fragmentation rules. This again 
shows that the relaxation time does not depend strongly 
on the absolute value of aoo- Instead, this indicates 



that the relaxation time is mainly determined by the 
collision rate between the particles. While a change in 
fragmentation rate, for example due to increased 7 dues 
not affect the relaxation time, an increased collision rate, 
due to an increased number of primary particles in the 
system seems to decrease the relaxation time significantly. 
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Figure 9: Histogram of the particle size distribution. Number of par- 
ticles versus the radius a for different values of the fractal dimension 
df oi the aggregates for (a) large-scale splitting (b) erosion splitting 
and (c) uniform splitting. 

Again, the size distributions (Fig. [T]) remain clearly dif- 
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ferent for the different fragmentation rules, independent of 
the total number of primary particles iVo in the system. 

3.5. Influence of the fractal dimension 

The last important system parameter that typically 
varies a lot in natural systems is the fractal dimension 
d/ of the aggregates 15]. To determine the influence of 
the fractal dimension of the aggregates, Uoo is computed 
for different values of d/. For all three fragmentation rules 
we find a drastic increase in the average number of pri- 
mary particles per aggregate in the steady state (see Fig. 
[8ja)) with increasing df. In the case of varying the fractal 
dimension this increase in ctoo is much more drastic than 
for the other parameters studied in the previous sections, 
ttoo increases by approximately a factor of 100 between 
df = 1.5 and df — 2.3. Initially, one might assume that 
the increase of the average number of primary particles 
per aggregate in the steady state is only due to the aggre- 
gates becoming more compact as the fractal dimension is 
increased and does not really reflect a change in the ag- 
gregate size. However, plotting the average radius of the 
aggregates in the steady state as a function df (inset in 
Fig. HJa)) shows that there is also a significant increase in 
the aggregate size with increasing df. 
This increase when varying the fractal dimension can be 
understood by looking at the stability condition for the ag- 
gregates (cf. Eq. ([T7| '). The stability curve defined by Eq. 
(jl7p becomes almost horizontal for larger aggregate sizes. 
The larger df the greater becomes the range of aggregate 
sizes where increasing the size has almost no effect on the 
stability (with the limit of d/ =3 where stability is in- 
dependent of the size) . Increasing the aggregate strength 
7 also leads to larger aggregates being stable at a given 
value of shear force. However, the increase in the range of 
stable aggregate sizes is much lower. 

We note that the relaxation time increases weakly with 
df for all splitting rules (see Fig. [HIb)). Increasing the 
fractal dimension leads to more compact aggregates, i.e. 
less overall volume occupied by aggregates and therefore 
smaller collision probabilities. This in turn increases the 
relaxation time. This effect appears strongest for erosion, 
where the relaxation time becomes very short for small 
values oi df whereas for large-scale fragmentation and uni- 
form fragmentation there seems to be a saturation of the 
relaxation time for values oi df < 2. 

Once again, the shapes of the particle size distributions 
retain their characteristic differences for the different frag- 
mentation rules. Aside from the increasing fluctuations in 
the distributions for increasing df, due to the decreasing 
number of aggregates in the system, the qualitative shape 
of the distribution remains characteristic for the fragmen- 
tation rule, independent of the value oi df. 
The range oidf that is observed in natural systems reaches 
even further than = 2.3, up to approximately 2.6. How- 
ever, this is not shown here, because due to the drastic 
increase in the average number of primary particles per 
aggregate only very few aggregates would remain for such 



high values oi df (if the other parameters remain fixed). 
Hence, no meaningful statistics or size distributions could 
be obtained. 

This is a computational limitation of the current model, 
due to the finite number of primary particles that can be 
studied. While it occurs for every parameter, e.g. increas- 
ing the aggregate strength 7 too far has the same effect, 
it is most pronounced for the fractal dimension. Small 
increases of the fractal dimension lead to proportionally 
much larger increases in the mean aggregate size, there- 
fore reducing the total number of aggregates in the system 
very quickly. 

4. Discussion and Conclusions 

To conclude this work, in the final section we provide 
a brief discussion of the limitations of the particle-based 
aggregation and fragmentation model proposed here and 
finish with a brief summary of our results. 

4-.1. Limitations of the Model 

In this work we showed the application of our 
particle-based model to the problem of aggregation and 
fragmentation of marine aggregates. We emphasize 
that the particle-based model introduced here is a very 
general model that can be applied to a wide range of 
other problems, for example in chemical engineering and 
has the potential to be a useful addition to the usual 
modelling approaches for aggregation and fragmentation. 
However, there are a number of limitations that should 
be pointed out. Some of these are due to physical aspects 
of the problem that are not yet fully understood and 
can therefore not be captured, others are mainly due to 
computational limitations. 

The main physical aspect of the problem that is still 
not fully understood are the details of the fragmentation 
mechanism. In particular for marine aggregates there still 
exists no satisfying microscopic theory for the fragmenta- 
tion process. The fragmentation model used in this work 
can therefore only be considered as a very simplified view 
on the problem and most likely does not capture many 
aspects of the real situation. 

Additionally, the equations of motion for fractal-like 
aggregates can be expected to be very different from the 
simple equations used here. The inclusion of the increased 
radius and effective density which was done in this work 
represents a very simple modification of the relevant 
forces, such as the drag forces acting on the particle and 
is unlikely to capture the full complexity of the problem. 
However, it is a useful approximation to consider the 
impact of a complex structure on the steady state size 
distribution. 

Furthermore, the description of marine aggregates as hav- 
ing a fractal-like structure is only a first approximation. 
Measurements have shown that when averaged over many 
different individual aggregates there exists a power-law 
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relationship between the aggregate size and mass. This 
is certainly not true for each individual aggregate, but 
has been found to be a reasonable model in many cases 

27| . In addition, the concept of a fractal dimension is 
usually only valid over a certain size range of aggregates. 
Aggregates that consist of only one or two primary 
particles generally do not have the same structure as a 
larger aggregate. In the context of this work it is possible 
to consider as primary particles the smallest fraction of an 
aggregate that still has the same power-law relationship 
as the large aggregates, i.e. as the fractal generator of the 
aggregate (see also the comment [s^). 
In addition, there are a number of aspects that are 
theoretically understood quite well, but can not be in- 
cluded in such a model due to computational limitations. 
First among these is the two- and three-way coupling 
between the particles and the surrounding fluid. Two-way 
coupling, i.e the feedback of the particles on the fluid can 
in principle be included but requires solving the equations 
of motion for the fluid together with the equations for 
the particles and therefore leads to a drastic increase in 
computational effort. In addition, as the particle radius 
r is assumed to be small, the feedback from the particle 
motion on the flow is usually not signiflcant [51] and 
can therefore be neglected. The inclusion of three-way 
particle-flow coupling, i.e. the interaction of particles 
through the fluid is usually a bigger challenge. In par- 
ticular at low Reynolds numbers this interaction affects 
the particle motion even at low particle concentrations 

52l |. While there exist a number of interaction models, 
that are able to compute this three-way coupling (see e.g. 
[H^) they are typically limited to systems with only a 
few particles because of the computational costs involved. 
While the modification of individual particle trajectories 
can be significant, the main effect for the collective 
aggregation and fragmentation dynamics of an ensemble 
of particle is typically a reduction of the collision rates. 
This can be approximated in a simple way by introducing 
a collision efficiency, i.e. a probability for aggregation 
upon collision. However, the introduction of such a 
collision efficiency does not qualitatively affect the results 
shown here. Additionally, the fluiddynamic interaction of 
permeable particles such as particles with a fractal-like 
structure discussed here are not understood very well. 
It is likely that flows through an aggregate can lead to 
very different interactions between fractal-like particles 
compared to solid spherical particles (see e.g [13, [s^), in 
particular such permeability effects may actually decrease 
the influence of hydrodynamic interactions between the 
aggregates compared to the case of completely solid 
particles fsil]. 

A further limitation for the applicability of the model is 
the number of primary particles that can be computa- 
tionally considered. The present computational capacities 
do not allow to apply this approach to large systems, e.g. 
models that are used to study aggregation, fragmentation 
and aggregate transport on spatial scales up to several 



hundred kilometers. For large systems a mean field 
approach is therefore much better suited. However, for 
many small systems and also for the principle study of 
the processes involved, this is not a severe limitation. 
Furthermore, as the model was tested with a simple 
3-d synthetic turbulent flow, it is an open question how 
representative the results are to draw general conclusions 
about the evolving steady state size distribution. In par- 
ticular intermittency effects and clustering of particles on 
the inertial scale of a real turbulent flow may significantly 
affect aggregation and fragmentation probabilities. In 
order to achieve more general statements it will therefore 
be necessary to study the model and the resulting size 
distribution for various, more realistic flows, for example 
using DNS simulations of real turbulent flows. Never- 
theless, the influence of different system parameters and 
fragmentation mechanisms has been tested and gives a 
detailed insight for this speciflc flow. 



4-2. Summary 

In the present study we described in detail a coupled 
model for advection, aggregation and fragmentation of in- 
dividual inertial particles with a fractal-like structure. We 
showed how typical properties of aggregation and frag- 
mentation processes can be incorporated. In particular, 
we introduced an approximate way, using modified ag- 
gregate sizes and effective densities, to account for the 
fractal-like structure that is common for aggregates in 
many natural systems. The model represents an alter- 
native approach to the mean field theory that is usually 
used to describe aggregation and fragmentation processes 
and was used to gain insights into principle behavior of 
fractal-like aggregates under different fragmentation mech- 
anisms. The model was parameterized for the case of a 
suspension of marine aggregates in the ocean, but can in 
principle be used in a wide range of applications such as 
cohesive sediment dynamics, the flocculation of biological 
cells or solid-liquid separation systems in chemical engi- 
neering [23,[ii,ti3|. 

We observed the development of a balance between ag- 
gregation and fragmentation, leading to a steady state. 
It was found that with increasing aggregate strength the 
mean aggregate size in steady state increases, whereas the 
relaxation time stays constant. With increasing fractal 
dimension the relaxation time towards steady state and 
the mean aggregate size in steady state increase. By con- 
trast, an increase in particle volume fraction decreases the 
relaxation time due to higher collision probabilities and 
increases the steady state mean size. In general, increased 
aggregation rates or decreased fragmentation rates lead to 
an increased mean aggregate size in steady state. The re- 
laxation time decreases for increasing aggregation rates, 
but does not change with decreasing fragmentation rates. 
In the context of our model different types of fragment size 
distributions can easily be tested and compared with each 
other. We compared numerical results for three commonly 



13 



used distributions of fragment sizes, large-scale fragmenta- 
tion where fragments typically have similar sizes, erosion, 
where one fragment is typically very small and uniform 
fragmentation, where all fragment sizes appear with the 
same probability. Large-scale fragmentation and erosion 
were treated as two separate processes to study the influ- 
ence of fragmentation at a distance ^ from the center of 
the aggregate on the size distributions in the steady state. 
In reality, aggregates will break with certain probabilities 
at a distance C from the center but there will not be two 
separate processes. Therefore, depending on the proba- 
bility distribution for fragmentation at a certain distance, 
one can expect different combinations of the steady state 
size distributions found in this work. 
One such combination, where fragmentation at all dis- 
tances ^ appears with the same probability was given by 
the uniform fragmentation rule. Uniform fragmentation 
leads to a distribution with a broad plateau for small ag- 
gregate sizes and an exponential tail towards larger ag- 
gregate sizes. A distribution with two different regimes 
evolves for erosion-like fragmentation. Large-scale frag- 
mentation leads to an exponential tail of the particle size 
distribution. Similar shapes of the size distribution of the 
aggregates for large-scale fragmentation have also been 
found in a number of theoretical and experimental stud- 
ies IB) 57 1, indicating that our model is able to reproduce 
the major features of such aggregation-fragmentation pro- 
cesses. Such an exponential tail has also been measured 
in field studies of marine aggregates 31| . This may indi- 
cate that large-scale fragmentation could be the primary 
mode of break-up for such aggregates, as has previously 
been discussed by e.g. Thomas et al. [18|. 
In all cases the steady state particle size distribution fol- 
lows a specific shape for each fragmentation rule. This in- 
dicates that the fragmentation process is most relevant for 
the shape of the distribution. The ratio of aggregation and 
fragmentation probabilities, mainly influenced by the ag- 
gregate strength, total particle volume fraction and fractal 
dimension, determines the mean aggregate size in steady 
state and the relaxation time. Out of these three parame- 
ters the fractal dimension has the strongest effect since it 
influences both aggregation and fragmentation probabili- 
ties. 

The influence of large-scale fragmentation versus erosion 
for marine aggregates has recently been studied numeri- 
cally and com par ed to experimental results in a work by 
Verney et al. [58[. They used a mean- field based Smolu- 
chowski equation approach and obtained results compara- 
ble to those of our model. Thus, as both model approaches 
lead to similar results the insight into fragmentation and 
fragment distributions provided by the perspective of our 
model can provide a useful addition to the understanding 
of aggregation and fragmentation processes. Additionally, 
in the particle-based model presented here particle inertia 
can be fully considered, while the correct incorporation 
of particle inertia into a mean field theory is still an un- 
solved problem. Hence, future model studies using this 



approach can lead to a better understanding of particle 
inertia effects in aggregation and fragmentation processes. 
Therefore, the model suggested here has the capability to 
be a powerful tool to investigate the validity of different 
approximative strategies in the formulation of a mean field 
theory. 
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